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ABSTRACT 

Studies of the Galactic population of radio pulsars have shown that their luminosity 
distribution appears to be lognormal in form. We investigate some of the consequences 
that occur when one applies this functional form to populations of pulsars in globular 
clusters. We use Bayesian methods to explore constraints on the mean and standard 
deviation of the luminosity function, as well as the total number of pulsars, given an 
observed sample of pulsars down to some limiting flux density, accounting for mea- 
surements of flux densities of individual pulsars as well as diffuse emission from the 
direction of the cluster. We apply our analysis to Terzan 5, 47 Tucanae and M 28, and 
demonstrate, under reasonable assumptions, that the number of potentially observ- 
able pulsars should be within 95.45% credible intervals of 133t5^\ 66^^^ and lOgtgf , 
respectively. Beaming considerations would increase the true population size by ap- 
proximately a factor of two. Future cluster pulsar discoveries would allow this method 
to be used to compare the luminosity function between different clusters. 

Key words: methods: numerical — methods: statistical — globular clusters: general 
— globular clusters: individual: Terzan 5, 47 Tucanae, M 28 — stars: neutron — 
pulsars: general 



Globular clusters are spherical collections of stars located 
throughout the halos of galaxies. Once thought to be com- 
posed entirely of old metal-poor population II stars, they 
are now also believed to precipitate during interactions or 
collisions of galaxies, therefor e containin g younger stars hav- 
ing higher metallicities (see IZepl l2003t ). The total masses 
of globular clusters ran ge up to the order of 10^ Mq (see 
iMevlan fc Heggielll997l ). and core stellar number densities 
reach 10^ pc"^. The high core densities lead to dynamical 
interactions between stellar systems that are found less com- 
monly in the Galactic plane. For example, globular clusters 
favour the formation of low-mass X-ray binaries (LMXBs) 
that ar e believed to be the progenitors of millisecond pulsars 
fMSPs: 'Alpa r et al.|[l983 ). and hence, the ratio of MSPs to 
normal pulsars is much larger in globular clusters than in 
the Galactic field. In addition, the MSPs in globular clus- 
ters tend to have higher eccentricities compared to their field 
counterparts, due to exchange or fly-by encounters. MSPs, 
due to their formation history, can be considered long-lived 



tracers of LMXBs, and therefore, constraints on the MSP 
content of globular clusters provide unique insights into bi- 
nary evolution and the integrated dynamical history of glob- 
ular clusters, while determining the radio luminosity func- 
tion of these pulsars helps shed light on the radio emission 
mechanism in action in these compact objects. 

Pulsar searches of globular cl usters have y ielded im- 
pressive returns in recent years (see lCamilo fc Ra sio 2005), 
with currently 143 pulsars known in 27 cluster^]. Of these, 
three clusters, Terzan 5, 47 Tucanae and M 28 are known 
to harbour more than 10 pulsars each, the most populous 
being Terzan 5 with 34 (Ransom S. M., private commu- 
nication). In this paper, we describe a Bayesian method 
that we have developed to compute an estimate of the true 
number of pulsars in a given cluster, given a set of ob- 
servations. There have been many attempts to constrain 
the population size of pulsars in all globular clusters in 
the G a laxy (see, for example, Ku lkarni, Narayan fc^ Romanil 
Il990l : iBagchi. Lorimer fc ChennamangalamI I201ll ). This 
work is different in that it treats clusters individually instead 
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of dealing with the total population. Bayesian approaches to 
constrain the pulsar population of individual globular clus- 
ters have been used specifically for the case of young (non- 
recycled) pulsars bv lBovles et alj (|201ll ). This work focuses 
on the entire radio pulsar content of the cluster - the ma- 
jority of which is made up of old (recycled) pulsars - and 
additionally, attempts to constrain luminosity function pa- 
rameters jointly with populati on size. 

Faucher-Giguere fc Kaspil (|2006 l have shown that the 
luminosity distribution of Galact ic pulsars appears to be 
lognormal in form. More recently, iBagchi et al.l |20l3) have 
verified that the observed luminosities of recycled pulsars 
in globular clusters are consistent with this result. Assum- 
ing, therefore, that there is no significant difference between 
the nature of Galactic and cluster populations, we investi- 
gate some of the consequences that occur when one applies 
this functional form to populations of pulsars in individual 
globular clusters. 

For a lognormal (base- 10) distribution of pulsar lumi- 
nosities, the luminosity function is given by 
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where L is the luminosity in mjy kpc'^, fj, is the mean and cr is 
the standard deviation of the distribution. We are interested 
in the situation where we observe n pulsars with luminosi- 
ties above some limiting luminosity Lmin. Given this sample 
of pulsars, we ask what constraints we can place on their lu- 
minosity function parameters, in addition to the potentially 
observable population size A*' (that is, the population of pul- 
sars beaming towards the Earth). Another way of thinking 
about this problem is that there is a family of luminosity 
function parameters and population sizes that is consistent 
with an observation of n pulsars above the luminosity limit 
of the survey, and we are analyzing the posterior probabil- 
ities of different members of this family given the observa- 
tions. 

This paper is organized as follows: In !j2l we describe 
our technique. In ^ we apply the technique to observa- 
tions of a few globular clusters to determine the constraints 
on the luminosity function parameters and population size. 
Later, we refine our results using a priori information on 
the luminosity function parameters to get a better estimate 
of the number of pulsars in those clusters. A summary and 
our conclusions are presented in 311 



2 BAYESIAN PARAMETER ESTIMATION 

Bayes' theorem (see I Wall fc Jenkins! [ioo^ : iGregorvl l2005h . 

for the purpose of parameter estimation, can be stated math- 
ematically as 



piO\D,M) 
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where is a set of parameters, D is some data and M 
is a model describing the parameters. In this notation, 
p{p\D, M) represents the probability of obtaining a set of pa- 
rameter values given the data and the model, and is termed 
the joint posterior probability density. Similarly, p{D\0, M) 
is the probability of having obtained the observed data, 
given the parameter values and the model, and is termed 



the likelihood, and p{0\M), the a priori probability dictated 
by the model, is termed the prior probability density. The 
denominator, p{D\M) is called the evidence, and is just a 
normalizing factor that can be dropped since we are only 
interested in relative probabilities, thereby giving 



pie\D, M) oc p{D\e, M)p{0\M). 



(3) 



In this paper, we use Bayes' theorem to find the joint pos- 
terior probability density functions of the model parameters 
^, a and N given some data. In our case, the data are the in- 
dividual pulsar flux densities that we call {Si}, the observed 
number of pulsars, n and the total diffuse flux density of the 
cluster, Sobs- 

Luminosity is a property intrinsic to pulsars, while flux 
density is the corresponding observable. The relationship be- 
tween the two quantities is given by 



47rr . 2,9. 
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where r is the distance to the pulsar, S is the pulse duty 
cycle, p is the radius of the pulsar emission cone, Smoan(j^) 
is the mean flux density of the pulsar as a function of ob- 
serving frequency and vi and 1/2 are the bounds of the 
frequency range over w hich the pulsar is observed (see 
iLorimer fc Krameill20 05lV Due to the uncertainty associated 
with the beam geometry, the values of S and p are not gen- 
erally reliable for luminosity calculations. Therefore, we use 
a simplified model of the luminosity, the 'pseudoluminosity', 
that is defined as = Sv r^ at a given frequency v (the 
1/ subscript on L and S will be dropped for the rest of the 
paper). As can be inferred from Equation Q and the afore- 
mentioned pseudoluminosity equation, the luminosity func- 
tion is inevitably corrupted by uncertainties in distance. To 
mitigate this, we decided to perform our analysis initially 
in terms of the measured fiux densities, and later, use a 
model of distance uncertainty to convert our results to the 
luminosity domain. We take the distance to all pulsars in a 
globular cluster to be the same. The lognormal in luminosity 
can then alternatively be written in terms of the flux den- 
sity. The probability of detecting a pulsar with flux density 
S in the range log S to log S + d(log S) is then given by a 
lognormal in S as 



p(log5) d(logS) = 



(log 
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where 5* is in mJy, and fj,s and as are the mean and standard 
deviation of the flux density distribution. The probability of 
observing a pulsar above the limit Smin is then 
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Our analysis involves computing three likelihoods in the 
flux domain based on three sets of data, computing the total 
likelihood as the product of these three likelihoods, convert- 
ing this fiux domain likelihood to the luminosity domain, 
and subsequently, applying priors to obtain the posterior. 
This procedure is depicted graphically in the block diagram 
of Figure [T] 
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Likelihood A: 
Pulsar flux densities 

n 
i=0 


LUieiihood B: 
Observed number of 
pulsars 
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Likeiihood C: 
Diffuse flux 







Totai likeliliood (flux): 

P(l0g Si, n, 5ob, I A', /(g, ffs. '^min) 
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Total likelihood (luminosity): 

p(log Sj, n, 5„bs I A', ff, o, 5^i„, r) 



Posterior: 
p(N, a, 5,„i„, r I log 5,, n, S„^,J 



Priors: 



p(N) 




P(SmJ 



Figure 1. Logical flow of the Bayesian analysis. The circle with the X sign symbolizes a multiplication operation. 



2.1 Using pulsar flux densities 

In the first step, we consider as data the measured fiux den- 
sities of the pulsars in the cluster under scrutiny, that we call 
{Si}. Ideally, the survey sensitivity limit Smin can be taken 
as another datum, but its exact value is not always known. 
The effects of radiometer noise, Doppler smearing, interfer- 
ence, and in some cases, interstellar scintillation, result in a 
distribution of Smin- We decided, therefore, to parametrize 
Smin. The likelihood of observing a set of pulsars with fluxes 
{Si} is represented as 

n 

PS ; 7 *5'min) , 

1=1 

where n is the number of observed pulsars in the cluster. 
Each term in this product is given by 

-j^ (log Sj-iis)^ 

Pi{logSi\ps,f^s,Smin) = y=e ""s (7) 

PobsO-sv27r 

where pobs is as given in Equation (|6]). This likelihood is 
represented as 'Likelihood A' in Figure [1] Uncertainties in 
the flux density measurements are not considered here, but 
it has to be borne in mind that it will have the effect of 
underestimating the credible intervals on our posteriors. 

2.2 Incorporating the number of observed pulsars 

To infer the total number of pulsars A'^ in the cluster, we fol- 
low Bovlos ct al. (2011) to take as likelihood the probability 
of observing n pulsars in a cluster with N pulsars, given by 
the binomial distribution 

p{n\N,^iS,(^S,Smin) = TT Fobs ( 1 ^ Pobs ) , (8) 

n\(l\ — n)! 



where Pobs is computed as in Equation This likelihood 
is shown as 'Likelihood B' in Figure [T] 

2.3 Considering diff^use emission 

Next, we incorporate information about the observed diffuse 
flux Sobs from the direction of the globular cluster. We as- 
sume that all radio emission is due to the pulsars in the clus- 
ter, including the both resolved pulsars and the unresolved 
background. Unlike standard pulsar searches, imaging the 
diffuse radio emission of a cluster to estimate the number 
of pulsars therein, is not affected by phenomena that cause 
pulse broadening, such as dispersion or scattering, or the 
fact that som e of the pulsars in the cluster are in acceler- 
ating frames l|Fruchter fc Gos3ll990l ). For the likelihood of 
measuring the diffuse flux Sobs, we choose 

p(Sobs|A^,Ms,<7s) = 3^ e "'^"d.ff" , (9) 

o"diffV 27r 

where Sdut is the expectation of the total diffuse flux of a 
cluster whose flux density distribution is a lognormal with 
parameters ps and as, and having N pulsars, and a^is is 
the standard deviation of this distribution. This likelihood is 
referred to as 'Likelihood C in Figure [T] Assuming that the 
cluster contains N pulsars, each having average luminosity, 

Sdiff = A(S), (10) 

and 

adiff = \/iV SD(S), (11) 
where the expectation of S is given by 

(S) = lO'^^+^'sWio)^ (^2) 
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and the standard deviation of S, 

SD(S) = lO^s + i's y^io-l Mio) _ i. (13) 

We do not consider the uncertainty in the diffuse flux mea- 
surement, and as mentioned in i\2.1l this has the effect of 
underestimating the credible intervals on our posteriors. 

The total likelihood, then is the product of the three 
likelihoods computed above, as shown below: 

p{l0gSi,n, SobslN, us, OS, Smin) 

71 

= n Pi (log •S'i I MS , 0"S , Smin) 



(14) 



X p{n\N, ^S, 0.5, Smin) 

X p(5'obs!A*', ^is,o"s). 



2.4 Transformation to luminosity domain 

The flux density distribution of pulsars in a globular cluster, 
although proportional to their luminosity distribution, is not 
suitable for comparing the populations in different clusters, 
as it depends on the distance to the cluster. It is, therefore, 
useful to transform the total likelihood obtained in the pre- 
vious subsection to the luminosity domain. We convert the 
total likelihood of Equation (|14|) to the luminosity domain in 
the following way. Since the pseudoluminosity equation can 
be written in terms of logarithms as log L = log S + 2 log r, 
where L is in mjy kpc^, S is in mjy, and r is in kpc, the 
means of the two distributions are related additively by the 
term 21ogr, while the standard deviations are the same. 
Taking into account the uncertainty in distance, we have a 
distribution of distances, p(r). The total likelihood in the 
luminosity domain is 



p(log Si,n, Sohs\N, fi, a, Smin, r) 

^p{logSi,n,SohB\N,fis,crs,Smin), 



(15) 



where fi and fis are related additively as mentioned above, 
and a and as are equal. The final joint posterior in lumi- 
nosity is then given by 



p{N,fi, a, Synin, r\ log Si,n, Sobn) 
cc p{logSi,n, SobslN, fi,a, Smin,r) 
X p{N) p(p) p(g) p{Sn,in) p{r). 



(16) 



The prior on A'' is taken to be uniform from n to 00, where 
the upper limit, for the sake of computation, would be an 
arbitrarily large value. We also use uniform priors on the 
model parameters /i and a. Due to the nature of the un- 
certainty in determining the exact value of Smin, we choose 
a uniform prior on it in the range (0, min(S'i)], where the 
upper limit is the flux density of the least bright pulsar in 
the cluster. The prior on r is taken to be a Gaussian. 

This joint posterior is integrated to obtain the following 
marginalized posteriors. For brevity, in these equations, D 



indicates the complete set of data used in this analysis. 
piN\D) = ffffp{N,fi,a,Sn.in,r\D) d^i da dSn^in dr 



(17a) 



3(/i|D) = jjjj p{N,fi,a,S,^in,r\D) dN da dSmin dr 

(17b) 



pia\D) = JJJJ P{N, M, o, Smin, r\D) dN dfi dS^in dr 

(17c) 

piS^in\D) = JJJJ p{N,fj,,a,Smin,r\D) dN dfi da dr 

(17d) 

p{r\D) = JJJJ A*, o, Smin, r\D) dN d^i da dSn^in 

(17e) 

The joint distribution of the parameters fi and a is given by 
marginalizing the joint posterior over A^, 5111111 and r, as 

p{li,a\D) = [[[ p{N,fi,a,Sn.in,r\D) dN dSn.in dr. (18) 



3 APPLICATIONS 

We applied our Bayesian techniqufl to the globular clusters 
Terzan 5, 47 Tucanae and M 28 (although the clusters we 
consider contain only recycled known pulsars, the analysis 
would remain the same even if there were young pulsars in 
the sample). The choice of clusters was based on the amount 
of data available. Terzan 5 is the cluster most-suited for this 
analysis due to the fact that it has a relatively large number 
of pulsars for which flux density measurements are avail- 
able. Although Terzan 5 has 34 known pulsars (Ransom S. 
M., private communication), we take n — 25, the number of 
pulsars for which we have flux density measurements. The 
flux densities of the individual puls ars were collected in a lit- 
erature survey (|Bagchi et al.ll20l"ll . and references therein), 
with the values relevant to this work reproduced in Table [T] 
For Terzan 5, the flux densities we used were s caled from 
those reporte d at 1 950 MHz by iRansom et"al] l|2005l 'l and 
iHessels et al] (|2006l ') to 1400 MHz using a s pectral index, 
a = —1.9, using the power la w S{iy) oc u". iHessels et al.l 
(|2007t) and iBagchi et al.l (|201ll ) discuss the choice of spec- 
tral index in detail. The observed diffuse flux density at 
1400 MHz is taken to be Sobs ~ 5.2 mJy (the sum of the 
diffuse flux and the fl uxes of point sources as reported by 
iFruchter fc Gossll2000l ). The priors used were formed in the 
following ways. The prior on A'' was chosen to be a uni- 
form distribution in [n, 500], w here the upper limit is 150% 
of the upper limit obtained by iBagchi et al.l (I2OII ) (using 
the values of fi and a as found bv lFaucher-Giguere fc Kaspil 
'2OO6) above their upper limit. We note that this prior is 
sufficiently wide to ensure that the posterior does not rail 
against the prior boundaries. We chose uniform distributions 
in the same range of fi and a as used bv lBagchl et al.1 (|201lD 

The software package that we developed to perform the anal- 
ysis described in this paper is available freely for download from 
|http://psrpop.ph ys . wvu . e du/gcbayes/| 
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as our priors. Survey sensitivity limits were not always avail- 
able, and additionally, due to a variety of factors mentioned 
in ii2.1l for all of our analyses, we took 5*111111 to be a uni- 
form distribution in the range (0, min(Si)]. The most recent 
measurement of the distance to Terzan 5, r — 5.5 ± 0.9 kpc 
l|Ortolani et al.ll2007l '), was used to model the distance prior. 
We modelled the distance prior as a Gaussian with mean 
5.5 kpc and standard deviation 0.9 kpc. Figure [2] shows 
the results of the analysis for Terzan 5. The mode of the 
marginalized posterior on TV, shown in Figure [IJb), is 35 
and the median with the surrounding 95.45% credible inter- 
val is 92lg4*. As can be seen from Figures [5{b), (Jl^c) and 
[2jd), the constraints we obtain on A*', jj, and a, respectively, 
are broad, due to a dearth of flux density measurements. 
The marginalized posterior on ^min, plotted in Figure [5Je), 
shows a strong preference for values away from and closer 
to that of the least bright pulsar observed. The main results 
are tabulated in Table [S] 

For 47 Tucanae and M 28, containing 14 and 9 pulsars 
each, the individual flux densities used are given in Table [1] 
We took gobs = 2.0 m jy (1400 MHz flux as reported by 
iMcConnell et al.ll2004h for 47 T ucanae, and gobs = 1 .8 mJy 
(1400 MHz flux as reported bv lKulkarni et al.lll990l ') for M 
28. The priors on N were taken to be uniform in the intervals 
[n, 225] for 47 Tucanae and [n, 400] for M 28, where the 
upper limits were computed in the same way as we did for 
Terzan 5. Priors on Smin were formed as in the case of Terzan 
5, in the range (0, m m(Si)] . Dist ance measurements from 
the Harris catalogue (|Harridll99d 2010 edition) were used 
to form the distance priors, with uncertainties taken to be 
15%. The main results for these clusters are tabulated in 
Table 

The value of can be further refined by considering 
possible dependences on other physical parameters of glob- 
ular clusters, and such a stu dy will appear in a future paper 
ijTurk et al.Jlin preparatimih . Note that A is the size of the 
population of pulsars in the cluster that are beaming to- 
wards the Earth. We can include the beaming fraction - 
the fraction of all pulsars beaming towards us - to refine 
this estimate. Uncertainties notwithstanding, the beaming 
fraction of millise cond pulsars is gen erally thought to be 
greater than 50% l|Kramer et al.|[T998l ). This, together with 
the fact that most pulsars in globular clusters are millisec- 
ond pulsars, imply that the true population size in a cluster 
is approximately a factor of two more than the potentially 
observable population size. 



3.1 Using prior information 

In the framework developed in the previous section, we use 
broad uniform (non-informative) priors for the mean and 
standard deviation of the lognormal. This lack of prior in- 
formation is apparent in Figure [S^b), where A' is not very 
well constrained. Prior infor mation can help bet ter constrain 
the parameters of inter est. iBovles et al.l ll201lll use models 
of Galactic pulsars from lRidlev fc Lorimeil l 2010l ) to narrow 
down /i to between —1.19 and —1.04, and a to the range 0.91 
to 0.98. We have chosen our priors on jj, and a to be uniform 
within these ranges. In this case, we have assumed that the 
luminosity function of globular clusters has the same form 



Table 1. Flux densities used in the analysis. 



Pulsar 


1400 MHz Flux Density (mJy) 


Terzan 5 


J1748-2446A 


1.91" 


J1748-2446C 


0.68" 


J1748-2446D 


0.08° 


J1748-2446E 


0.09° 


J1748-2446F 


0.07" 


J1748-2446G 


0.03" 


J1748-2446H 


0.03" 


J1748-2446I 


0.05" 


J1748-2446J 


0.04" 


J1748— 2446K 


0.08" 


J1748-2446L 


0.08" 


J1748-2446M 


0.06" 


J1748-2446N 


0.10" 


J1748-24460 


0.23" 


J1748-2446P 


0.14" 


J1748-2446Q 


0.05" 


J1748-2446R 


0.02" 


J1748-2446S 


0.03" 


J1748-2446T 


0.04" 


J1748-2446U 


0.03" 


J1748-2446V 


0.13" 


J1748-2446W 


0.04" 


J1748-2446X 


0.03" 


J1748-2446Y 


0.03" 


J1748-2446ad 


0.15'' 


47 Tucanae 


J0023-7204C 


0.36'= 


J0024-7204D 


0.22= 


J0024-7205E 


0.21= 


J0024-7204F 


0.15= 


J0024-7204G 


0.05= 


J0024-7204H 


0.09= 


J0024- 72041 


0.09= 


J0023-7203J 


0.54= 


J0024-7204L 


0.04= 


J0023-7205M 


0.07= 


J0024-7204N 


0.03= 


J0024-7204O 


0.10= 


J0024-7204Q 


0.05= 


J0024-7203U 


0.06= 


M 28 


B1821-24A 


0.94"* 


J1824-2452B 


0.07'' 


J1824-2452C 


0.17'^ 


J1824-2452D 


0.05'^ 


J1824-2452E 


0.06'^ 


J1824-2452F 


0.08'^ 


J1824-2452G 


0.05"^ 


J1824-2452H 


0.06'^ 


J1824-2452J 


0.07'* 



Based on values reported by I Ransom et aL I hood ) 
^ Based on values rep orted bv lHessels et aL I ll2006l) 
ICamilo et al] l l200oi') 
iBeeinI 1I2OO6I ') 



as those in the Galactic field. This a ssumption appears t o 
be justified by reasons elaborated in iBagchi et al.l l|201lf ). 
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Figure 2. Results of the Bayesian analysis for Terzan 5, with {Si} as give n in Table [T] n = 2 5, and Sohs = 5-2 mjy. This analysis 
was run with wide priors on ^ and a, with the ranges equal to those used bv lBagchi et alj l|201lh (their Figure 2). (a) depicts the joint 
posterior on fi and cr, marginalized over N, S^nin a^nd ^ given in Equation II18I I: (b) is the marginalized posterior for as given in 
Equation l|17fe .'). with a mode of 35 and a median of 92. The x-axis is plotted in log scale for clarity; (c) is the marginalized posterior for fi 
as given in Equation Ill7b ) with a mode of —0.35 and a median of —0.9; (d) is the marginalized posterior for a as given in Equation I ll7b ) 
with a mode of 1.0 and a median of 0.9; (e) is the marginalized posterior for 5niin as given in Equation ((iTjl) with both mode and median 
equal to 0.02 mJy; (f) is the marginalized posterior for r as given in Equation l UTfe ). with a mode of 5.56 kpc and a median of 5.45 kpc. 
The hatching indicates regions that lie outside a 95.45% credible interval. 



Applying the Bayesian analysis over this narrower range of 
jj, and (T for Terzan 5 results in much tighter constraints on 
A'^ as seen in Figure [31 a), in which the mode of the distribu- 
tion is 123 and the median and a 95.45% credible interval 



is 133: 



The analysis was also performed for 47 Tucanae 



and M 28, the results of which are given in FiguresO^b) and 
[3jc), respectively. For 47 Tucanae, the mode of is 61 and 
the median with the surrounding 95.45% credible interval is 
66I3J. For M 28, the mode is 94 and the median with credi- 
ble interval is 109lfi^^. Our resu lt for Terzan 5 is consistent 
with that of*Bag chi et al.l (l201lh . In the case of 47 Tucanae 
and M 28, there is partial, yet considerable overlap between 
our credible ran ges and the c o rrespo nding confidence inter- 
vals quoted by iBagchi et al.l (I2011I). For 47 Tucanae , our 
results agree well with those of iGrindlav et al.| (|2002l '). i.e., 
35-90 MSPs with X-ray luminosities above 10^° ergs s~^. 
The results of our analyses are tabulated in Table [2l 



Table 2. Median and credible intervals from the various analyses 
presented in this paper. We note here that, in addition to the 
sources of error mentioned in the text, the values of fi and a 
presented here are also affected by the fact that computations 
are discrete and hence use a finite number of steps. 
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Figure 3. Posteriors on after applying tlie iBovl es et al ] l|201ll 'l priors on fj, and cr: (a) Terzan 5, with n = 25 and 5obs = 5.2 mjy. 
Tile median witii tiie surrounding 95.45% credible interval of N is 133^gg^; (b) 47 Tucanae, with n = 14 and 5obs = 2.0 mJy. The 
median with credible interval is 66^3^; (c) M 28, with n = 9 and Saha = 1-8 mJy. The median with credible interval is 109^g3^. The flux 
densities of the individual pulsars, {Si}, used in this analysis are given in TablefT] 



4 SUMMARY AND CONCLUSIONS 



We have developed a Bayesian technique to constrain the lu- 
minosity function parameters and population size of pulsars 
in individual globular clusters, given a data set that con- 
sists of the number of observed pulsars, the flux densities of 
the individual pulsars in the duster and the total diffuse flux 
emission from the direction of the globular cluster, assuming 
a lognormal luminosity function. We have applied our anal- 
ysis to a few globular clusters and have demonstrated the 
utility of this technique in constraining the aforementioned 
parameters. 

Our technique is applied in two different ways - first, 
with no prior information, and second, assuming prior 
knowledge of the possible range of /i and a. As shown for 
Terzan 5, the results for the first approach do not constrain 
N, fi OT a very well, but the lat ter two do exhibit con sis- 
tency with the values found by iFauc hcr-Gigucrc & Kaspil 
l|2006l ') and lBagchi et al.l (|201ll 'l. For the second approach in 
which we assume prior information to bound jU and a, the 
priors help better constrain the total number of pulsars in 
the cluster. 

The technique we have developed here should prove use- 
ful in future studies of the globular cluster luminosity func- 
tion where ongoing and future pulsar surveys are expected 
to provide a substantial increase in the observed populations 
of pulsars in many clusters. In particular, we anticipate that 
the increased amount of data would enable us to constrain 
the distributions of /i and a independently (i.e. without the 
need to assume prior information from the Galactic pulsar 
population). Further interferometric measurements of the 
diffuse radio flux in many globular clusters could provide 
improved constraints on fi and a by measuring the flux con- 
tribution from the individually unresolvable population of 
pulsars. 
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